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Abstract. Snake "resonances" are classified in terms of the invariant spin field and the amplitude 
dependent spin tune. Exactly at snake "resonance" there is no continuous invariant spin field at most 
orbital amplitudes. 

PROLOGUE 

This is an extended version of the paper with the same title published in the proceedings 
of the conference SPIN2002 [1]. A key aspect of the original paper was that the invariant 
spin field for snake "resonances" is irreducibly discontinuous at most orbital amplitudes. 
However, details were omitted owing to the page limit. In the meantime other papers 
[2, 3] have appeared which discuss the invariant spin field at snake "resonances" and it 
has become clear that it would be useful to extend [1] to give more details. 

We begin by presenting a slightly polished version of the original paper. Then the 
additional material is presented as an addendum. The citations are also updated. 



INTRODUCTION 

Spin motion in storage rings and circular accelerators is most elegantly systematised in 
terms of the invariant spin field (ISF) and the amplitude dependent spin tune (ADST). 
Here we apply them in the context of snake "resonances". We begin by briefly recapitu- 
lating some necessary basic ideas. For more details see [4, 5, 6, 7]. 

Spin motion in electric and magnetic fields at the 6-dimensional phase space point z 
and position s around the ring, is described by the T-BMT precession equation dS/ds — 
€l(j£;s) x S [8, 4] where S is the spin expectation value ("the spin") in the rest frame of 
the particle and €l(j£;s) contains the electric and magnetic fields in the laboratory. The 
ISF, denoted by n(z;s), is a 3-vector field of unit length obeying the T-BMT equation 
along particle orbits (z(s);s) and fulfilling the periodicity condition n(z;s + C) = n(z;s) 
where C is the circumference. Thus n(M(z;s); s + C) = n(M(z;s);s) = R 3x3 (z;s)n(z;s) 
where M(z; s) is the new phase space vector after one turn starting at z and s and R 3x3 (z; s) 

is the corresponding spin transfer matrix. The scalar product 7 S = S ■ h is invariant along 
an orbit, since both vectors obey the T-BMT equation. Thus with respect to the local h 
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the motion of S is simply a precession around n. The field n can be constructed at each 
reference energy where it exists without reference to individual spins. 

The chief aspects of the ISF are that: 1) For a turn-to-turn invariant particle distri- 
bution in phase space, a distribution of spins initially aligned along the ISF remains 
invariant (in equilibrium) from turn-to-turn, 2) for integrable orbital motion and away 
from orbital resonances the ISF determines the maximum attainable time averaged po- 
larisation P Um = | (n(z; s) ) | on a phase space torus at each s, where () denotes the average 
over the orbital phases, 3) under appropriate conditions J s is an adiabatic invariant while 
system parameters such as the reference energy are slowly varied, 4) it provides the 
main axis for orthonormal coordinate systems at each point in phase space which serve 
to define the ADST which in turn is used to define the concept of spin-orbit resonance. 

These coordinate systems are constructed by attaching two other unit vectors u\(z;s) 
and U2(z;s) to all such that the sets (u\,U2,n) are orthonormal. Like n, the fields u\ 
and U2 are 1 -turn periodic in s: Ui(z;s + C) = Ui(z;s) for i'e{l,2}. With the basis vectors 
u\ and U2 we can quantify the rate of the above mentioned spin precession around ft: 

— » 

it is the rate of rotation of the projection of S onto the u\, U2 plane. Except on or close 
to orbital resonance, the fields ui(z;s) and foil', s) can be chosen so that the rate of 
precession is constant and independent of the orbital phases [9, 5, 6, 7]. The number of 
precessions per turn "measured" in this way is called the spin tune. The spin tune, V S (J), 
depends only on the orbital amplitudes (actions) J, hence the name ADST. The choice 
of some ui(z;s) and &2(z;s) satisfying the condition u~i(z;s + C) = Ui(z;s) for /e{l,2} 
is not unique. An infinity of others can be chosen by suitable rotations of the k,- around 
ft. These lead to the equivalence class of spin tunes obtained by the transformation: 
V s (7) =>■ V S (J) + Iq + l\Qi + I2Q2 + I3Q3 for any integers / where the Q(J)'s are the 
tunes on a torus of integrable orbital motion l . The ADST provides a way to quantify 
the degree of coherence between the spin and orbital motion and thereby predict how 
strongly the electric and magnetic fields along particle orbits disturb spins. In particular, 
the spin motion can become very erratic close to the spin-orbit resonance condition 
V S (J) = mo + miQi+ m2<22 + ^3<23 where the m's are integers. At these resonances 
the ISF can spread out so that P lim is very small. Examples of the behaviour of P Um near 
spin-orbit resonance and the application of a generalised Froissart-Stora description of 
the breaking of the adiabatic invariance of 7 S while crossing resonances during variation 
of system parameters can be found in [6, 7, 10, 1 1]. Note that: 1) the resonance condition 
is not expressed in terms of the spin tune v(0) on the closed orbit, 2) a "tune" describing 
spin motion but depending on orbital phases could not be meaningful in the spin-orbit 
resonance condition, 3) if the system is on spin-orbit resonance for one spin tune of 
the equivalence class, it is on resonance for all others. In general u\ and U2 do not obey 
the T-BMT equation along an orbit (z(s);s). But at spin-orbit resonance, they can be 
chosen so that a spin S is at rest in its local (ui,U2,n) system. Then ui(z;s) and U2(z',s) 
do obey the T-BMT equation so that the ISF n(z;s) is not unique. 

Nowadays we emphasise the utility of the ISF for defining equilibrium spin distri- 
butions. However, it was originally introduced for bringing the combined semiclassi- 



For a recent detailed discussion of these concepts, see [3]. 
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cal Hamiltonian of spin-orbit motion into action-angle form for calculating the effects 
of synchrotron radiation [12]. The initial Hamiltonian is written as H s = ^{Q\J\ + 

22-^2 + 23/3) + & • S. By viewing the spin motion in the (u\ , u% n) systems, a new Hamil- 
tonian in full action-angle form Hf_ Q = ^(Q[J[ + + Q3J3) + ^V s (f)J s is ob- 
tained which is valid at first order in ft [9]. This emphasises again that, as with all action- 
angle formulations, the spin frequency cannot depend on orbital phases. Moreover, it is 
easy to show that at orbital resonance, ( i.e. ko + k\ Q\ + &2<22 + ^3<23 =0 for suitable 
integers k) the "diagonalisation" of the Hamiltonian (i.e. finding the «i,«2) might not 
be possible [9]. Thus at orbital resonance the ADST may not exist. On the other hand, 
one avoids running a machine on such resonances. The spin tune on the closed orbit 
Vo = v s (0) always exists and so does ho{s) = n(0;s). 

For our present purposes there are two kinds of orbital resonances: resonances where 
at least one of the Q's is irrational and those where all are rational. We write the rational 
tunes as Qi = cii/bi (i = 1,2,3) where the a ; and bi are integers. Then for the second 
type, the orbit is periodic over c turns where c is the lowest common multiple of the bi. 
This opens the possibility that in this case the ISF at each (z, s) can be obtained (up to a 
sign) as the unit length real eigenvector of the c-turn spin map (c.f. the calculation of ho 
from the 1-turn spin map on the closed orbit.). However, the corresponding eigentune 
cv c extracted from the complex eigenvalues A c = e ±2mcVc , depends in general on the 
orbital phases at the starting z. Thus in general v c is not a spin tune and should not be 
so named [13]. Nevertheless if c is very large the dependence of V c on the phases can be 
very weak so that it can approximate well the ADST of nearby irrational tunes. This is 
expected heuristically since the influence of the starting phase can be diluted on forming 
the spin map for a large number of turns. At non-zero amplitudes, both for irrational or 
rational <2's, the eigentune of the 1-turn spin map usually has no physical significance. 
Of course, it normally depends on the orbital phases and the corresponding eigenvectors 
are normally not even solutions of the T-BMT equation. 

For non-resonant orbital tunes, the spin tune can be obtained using the SODOM-II 
algorithm [14] whereby spin motion is written in terms of two component spinors and 
SU(2) spin transfer matrices. The functional equation n(M(z;s);s) = R 3xi (z;s)h(z;s) 
is then expressed in terms of a Fourier representation, w.r.t. the orbital phases, of the 
spinors and of the 1-turn SU(2) matrices. The spin tune appears as the set of eigentunes 
of an eigen problem for Fourier components and h is reconstructed from the Fourier 
eigenvectors. SODOM-II delivers the whole spin field on the torus J at the chosen s. 



THE SINGLE RESONANCE MODEL 

In perfectly aligned flat rings with no solenoids, no is vertical and Vo is in the equiv- 
alence class containing ayo where 70 is the Lorentz factor on the closed orbit and 
a is the gyromagnetic anomaly of the particle. In the absence of skew quadrupoles, 
the primary disturbance to spin is then from the radial magnetic fields along verti- 
cal betatron trajectories. The disturbance can be very strong and the polarisation can 
fall if the particles are accelerated through the condition ajo = K = ko ± Qi where 
mode 2 is vertical motion. This can be understood in terms of the "single resonance 
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model" (SRM) whereby a rotating wave approximation is made in which the contri- 
bution to Q. from the radial fields along the orbit is dominated by the Fourier har- 
monic at K with strength efa). The SRM can be solved exactly and the ISF is given 
by [15] n($2) = sgn(<5) (8e2 + e(eicos02 + e^sin^)) /V8 2 + e 2 where 5 — ajo — k, 
§2 is the orbital phase, (ei,^,^) are horizontal, vertical and longitudinal unit vec- 
tors and the convention h ■ £2 > is used. The tilt of h away from the vertical ho is 
I arcsin(e/v / S 2 + e 2 ) | so that it is 90° at 8 = for non-zero £. At large the equilib- 
rium polarisation directions h(J2, fas), are almost parallel to fio(s) but during acceler- 
ation through 6 = 0, h varies strongly and the polarisation will change if the adiabatic 
invariance of J s violated. The change in J s for acceleration through 8 = is given by 
the Froissart-Stora formula. The ADST which reduces to afo on the closed orbit is 
v s = sgn(<5)V <5 2 + e 2 + k. Note that the condition 8 = is not the spin-orbit resonance 
condition. On the contrary, as 8 passes through zero v s jumps by 2e with our convention 
for h and avoids fulfilling the true resonance condition: for particles with non-zero £, 
afo is just a parameter. In this simple model v s exists and is well defined near spin- 
orbit resonances for all Q2. This is also true in more general cases if orbital resonance is 
avoided. 



THE SINGLE RESONANCE MODEL WITH A PAIR OF SIBERIAN 

SNAKES 

Snake "resonances" 

Polarisation loss while accelerating through 5 = can be reduced by installing pairs 
of Siberian Snakes, magnet systems which rotate spins by n independently of z around 
a "snake axis" in the machine plane. For example, one puts two snakes at diametrically 
opposite points on the ring. Then ho ■ £2 = +1 in one half ring and —1 in the other. With 
the snake axes relatively at 90°, Vo is in the equivalence class containing 1/2 for all Jo- 
For calculations one often represents the snakes as elements of zero length ("pointlike 
snakes"). Then if, in addition, the effect of vertical betatron motion is described by the 
SRM, and orbital resonances are avoided, calculations with SODOM-II, perturbation 
theory [16] and the treatment in [13] suggest that V s (/2) is in the equivalence class 
containing 1/2 too, independently of 70 but also of J2. Thus for Q2 sufficiently away 
from 1/2 no spin-orbit resonances V S (J2) = ^0 ± Qi are crossed during acceleration 
through 5 = and the polarisation can be preserved. This is confirmed by tracking 
calculations. However, such calculations and analytical work show that the polarisation 
can still be lost if the fractional part of Q2, [Qi], is ^2/2^2 where here, and later, 02 and 
h>2 are odd positive integers with 0,2 < 2t>2 and where here and later the brackets [...] 
are used to signal the fractional part of a number. This is the so-called "snake resonance 
phenomenon" and it also has practical consequences [17, 18, 19], especially for small l>2- 
Such a [Q2] fits the condition 1/2 = (1 — a2)/2 + &2[<22]- But calculations (see below) 
show that exactly at [Q2] = ^2/2^2 the ADST may not exist. If it doesn't, it isn't in 
the equivalence class for 1/2. Then we are not dealing with a conventional resonance 
V s (^2) = (1 — «2)/2 + &2[<22] and the term resonance is inappropriate. Depolarisation 
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in this model has also been attributed to the fact that for non-zero J2 the eigentune of 
the 1-turn spin map, which depends on (j) 2 , is 1/2 at some values of 2 [17]. However, 
such a quantity does not describe spin-orbit coherence. Snake "resonances" are usually 
associated with acceleration but it has been helpful in other circumstances [6, 7, 10] to 
begin by studying the static properties of the system, namely with the ISF. We now do 
that for the SRM with two snakes for representative, parameters. 



Numerical study 

Figure 1 shows P lim (just before a snake) and v s for 25000 equally spaced [Q2]'s 
between and 0.5 for e = 0.4 and 8 = 0. At each [Q2], n is calculated by stroboscopic 
averaging [4] (< 25 10 6 turns) at 500 equally spaced 02 in the range — 2n and P lim 
is obtained by averaging over these 02. The ADST is obtained from SODOM-II. If the 
ADST exists SODOM-II delivers a part of the equivalence class, namely the spectrum 
[±0.5 ± I2Q2] for a range of contiguous even I2 restricted by the necessarily finite size 
of the matrix of Fourier coefficients. Only even I2 are allowed by the algorithm. For 
irrational Q2 the range of I2 is large. For rational Q2 the spectrum can include ±0.5 
but is otherwise highly degenerate or contains none or just a very few of the required 
members [±0.5 ± hQi]- Thus the existence of an ADST is easily checked. The central 
horizontal row of points in figure 1 shows the common member ±0.5 of the equivalence 
class of the ADST at the values of [Q2] where the ADST exists. There is an ADST 
for most [Q2] 's used. The first row of dots up from the bottom marks [Q2] values 
where there is no ADST. As expected, these are all at rational [Q2]'s such as 1/5, 1/4, 
2/5 ... or ^2/2^2 = 1/6,3/14,3/10... and the [cv c ] computed for these [Q2] show 02 
dependence. The curved line shows P lim and the second row of dots from the bottom 
marks [Q2] values where the ISF obtained by stroboscopic averaging did not converge 
for all phases. These coincide with sharp dips in P lim and are at or near [Q2] = #2/2^2, 
i.e in the snake "resonance" subset of the [22] ' s m the first row. Thus snake "resonance" 
is already a static phenomenon. Near such [Cfel's, the ISF, which for just one orbital 
mode is a closed curve in three dimensions, becomes extremely complicated as h strives 
to satisfy its defining conditions. Right at [Q2] = ^2/2^2 the nonconvergence occurs 
at [02/2tt] = j/2^2 for integers j = 1,...,2^2 and, moreover, the ISF is discontinuous 
at these phases 2 . For [Q2] = ^2/4^2 («2 < 4Z?2)> P lim an d the ISF show no special 
behaviour. These observations are consistent with the perturbative result [7] that for 
mid-plane symmetric systems, h should be well behaved near even n%2 but may show 
exotic behaviour close to odd n%2 = i>2. As expected, P lim and the ISF also show no 
special behaviour for [Q2] =^2/^2 («2 < ^2)- Some snake "resonances" such as that at 
[Q2] = 1/30 are narrower than 0.00002 in [Q2] and are missed in this scan. P Um also has 
several dips at values of [Q2] (e.g. at 0.341) which appear to have no special significance, 
but which should still be avoided at storage. The results for 0.5 < [Q2] < 1.0 are the 
reflection in 0.5 of the curves and points shown. Qualitatively similar results are obtained 



Note that in [3] it was convenient to require that the magnetic and electric fields and the ISF were smooth 
in (j) and s. Here we drop that requirement since we are dealing with models with pointlike snakes. 
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with equally distributed odd pairs of snakes set to give Vo = 1/2. The ISF and P lim usually 
vary significantly with s. 
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FIGURE 1. P lim (left axis) and a component of the ADST (right axis) for the SRM with 8 = 0, £ = 0.4 
and with 2 Siberian Snakes with axes at 90° and 0°. 



0.5 



E 





Radial component 
Vertical component 
Longitudinal component 








AAA 


AA 













0.5 1 

FIGURE 2. The three components of n(02) for the SRM with 2 Siberian Snakes with axes at 90° and 
0° and for [Q2] — 1/6. Viewing point: just before a snake. 5 = and e = 0.4. 
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Addendum 



A typical example for an ISF at a snake "resonance" is shown in figure 2. This shows 
the components of the ISF for [Q2] = 1/6 and £ = 0.4. In this case they are obtained by 
using the 6-turn spin map to calculate the vector h in the range < [02/27T] < 1/2&2, 
namely to 60 degrees, while applying the constraint that h should be continuous in 
02. Then the h for each 02 in this range is transported with the 1-turn spin map for 
six, or more, turns. One sees that h changes sign at values of 02 which are multiples of 
60 degrees so that the ISF is discontinuous as advertised. The six sets of stray points 
at [02/2tt] = 1/12, 1/4. . . are at phases 02 where the 6-turn spin map is the identity. 
At these phases h obtained in this way is arbitrary and the algorithm delivers values 
dominated by numerical noise. Of course, if the ISF is represented as the locus of points 
on the surface of the unit 2-sphere, this ISF gives disjoint segments. It is also easy 
to demonstrate that the positions of the discontinuities can be shifted. Thus at snake 
"resonances" the invariant spin field is not only discontinuous but also non-unique and 
to an extent which goes far beyond the non-uniqueness at the phases 02 at which the 
6-turn spin map is the identity. The ISF shows analogous behaviour at other snake 
"resonances". This is in stark contrast to the case, say, of the pure SRM. Note that 
discontinuities would not be allowed for irrational [Q2] , but that they are not prohibited 
for rational [Q2] • Moreover, theorems on uniqueness of the ISF require the existence of 
a spin tune and that the system is away from orbital resonance. Of course, for snake 
"resonances" with very high values of h>2, the ISF and the corresponding configuration 
of equilibrium polarisation is very complicated. It is then far from clear whether the 
ISF is a useful concept for these simple models involving just one plane of orbital 
motion and singular, i.e., non-physical fields, although, of course, such models have 
been instrumental in presaging the loss of polarisation observed in real storage rings 
[17,18,19]. 

The ISF obtained from the real eigenvector of the multi-turn spin map is also non- 
unique for [Q2] = ^2/4^2, namely at those [02/2tt] where the 4^2-turn spin map is the 
identity. Moreover, additional discontinuities in the sign can be added by hand, thereby 
enhancing the choice of ISF's. But in contrast to the case of snake "resonances" such 
discontinuities remain optional. 

Further aspects of these matters will be reported elsewhere. 

It is instructive to compare the discontinuous curves in figure 2 with the smooth curves 
in figures 7 and 8 in [2] 3 . These are also said to represent ISF's at snake "resonances". 
The vectors corresponding to the curves in figure 2 satisfy the T-BMT equation by 
construction and they are single valued in [02/2tt] as required for an ISF. However, 
if the vectors for the curves in figures 7 and 8 in [2] are transported according to the 
T-BMT equation, they are not single valued in [02/2^]. Alternatively, if those curves 
are taken to represent single-valued functions of [02/2 ft], as depicted, then they do not 
represent spin motion, i.e., motion according to the T-BMT equation. Either way, the 
curves in those figures do not represent ISF's at snake "resonances". 



It is also shown in [2] that the curve in figure 1 can be reproduced by the MILES algorithm by the same 
author. MILES, like SODOM-II, is based on SU(2) and Fourier expansions. 
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SUMMARY 



A snake "resonance" is at root a static phenomenon characterised by an invariant spin 
field which, for the simple models discussed here, is irreducibly discontinuous in 02 
for most orbital amplitudes. Moreover, on and near snake "resonance", there is no 
amplitude dependent spin tune so that the snake "resonances" of these models are 
not simple spin-orbit resonances. The mechanism, in terms of J s , for polarisation loss 
during acceleration through 8 = at and near such [Q2]'s is under study. 

We thank K. Heinemann, G. H. Hoffstaetter and J.A. Ellison for useful discussions. 
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